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Abstract: We present an up-to-date global analysis of solar, atmospheric, reactor and 
accelerator neutrino data in the framework of three-neutrino oscillations. We discuss in 
detail the statistical significance of the observed "hint" of non-zero #13 in the solar sector at 
the light of the latest experimental advances, such as the Borexino spectral data, the lower 
value of Gallium rate recently measured in SAGE, and the low energy threshold analysis of 
the combined SNO phase I and phase II. We also study the robustness of the results under 
changes of the inputs such as the choice of solar model fluxes and a possible modification 
of the Gallium capture cross-section as proposed by SAGE. In the atmospheric sector 
we focus on the latest results for v e appearance from MINOS and on the recent Super- 
Kamiokande results from the combined phases I, II and III, and we discuss their impact 
on the determination of #13. Finally, we combine all the data into a global analysis and 
determine the presently allowed ranges of masses and mixing. 
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1. Introduction 

It is now an established fact that neutrinos are massive and leptonic flavors are not sym- 
metries of Nature [1,2]. In the last decade this picture has become fully proved thanks to 
the upcoming of a set of precise experiments. In particular, the results obtained with solar 
and atmospheric neutrinos have been confirmed in experiments using terrestrial beams: 
neutrinos produced in nuclear reactors and accelerators facilities have been detected at 
distances of the order of hundreds of kilometers [3] . 

The minimum joint description of all the neutrino data requires mixing among all the 
three known neutrinos (v e , v„, v T ), which can be expressed as quantum superpositions of 
three massive states (i = 1, 2, 3) with masses mj. This implies the presence of a leptonic 
mixing matrix in the weak charged current interactions [4, 5] which can be parametrized 
as: 

C12C13 S12C13 si 3 e~ lScp \ 

U = I -S12C23 - ci2Si 3 s 2 3e wScp Ci 2 C 23 - Sl2Sl3S23e i<5cP C13S23 



(1.1) 



■S12S23 - C 12 S 13 C2 3 e' l5cp -C12S23 - Sl2Sl3C23e i<5cP C13C23 ) 



where Cij = cos Oij and Sij = sin^. In addition to the Dirac-type phase Sep, analogous 
to that of the quark sector, there are two physical phases associated to the Majorana 
character of neutrinos, which however are not relevant for neutrino oscillations [6, 7] and 
are therefore omitted on the present work. Given the observed hierarchy between the solar 
and atmospheric mass-squared splittings there are two possible non-equivalent orderings 
for the mass eigenvalues, which are conventionally chosen as 

Am^ < (Am^ ~ Am^) with (Am| x > 0) ; (1.2) 
Amli < lAm^ ~ Amg 2 | with (Am^ < 0) . (1.3) 
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As it is customary we refer to the first option, Eq. (1.2), as the normal scheme, and to 
the second one, Eq. (1.3), as the inverted scheme; in this form they correspond to the two 
possible choices of the sign of Am^. In this convention the angles 9{j can be taken without 
loss of generality to lie in the first quadrant, 0y € [0, 7r/2], and the CP phase <5cp £ [0, 2tt]. 

Within this context, Am^, (Am^l, #12, and #23 are relatively well determined [3,8- 
10], while only an upper bound is derived for the mixing angle #13 and barely nothing is 
known on the CP phase Sqp and on the sign of Am^. Apart from the importance from a 
theoretical point of view, establishing whether #13 is zero or not is an essential step in the 
development of the search strategies for the upcoming experiments. Once we know that the 
two mixing angles 612 and #23 are relatively large, the possibility of experimentally accessing 
leptonic CP violation crucially depends on the value of the angle #13. Also a non-zero #13 is a 
fundamental ingredient for a feasible determination of the neutrino mass ordering. For this 
reason, it is a main objective of upcoming reactor and accelerator experiments to directly 
measure this parameter. In this respect, Refs. [11,12] pointed out that two independent 
hints in favor of a non-zero value of 6*13 emerge from the combination of solar and long- 
baseline (LBL) reactor data as well as from the combination of atmospheric, short-baseline 
reactor and LBL accelerator data. Since these signals are the results of synergies between 
different data samples, it is particularly important to verify their stability with respect to 
new experimental data as well as to variations on the assumptions in the analysis. 

In this work we present the results of an up-to-date global analysis of solar, atmo- 
spheric, reactor and LBL accelerator neutrino data in the context of three-neutrino oscil- 
lations. In Sec. 2 we focus on the solar sector and we assess the stability of the oscillation 
parameters with respect to the inclusion of the new experimental results and theoretical 
advances which have become public during the last year. In this context, we find that 
many of these changes indeed lower the statistical significance of a non-zero value of #13. 
In Sec. 3 we do the same in the atmospheric sector, with particular emphasis on the recent 
v e appearance results from MINOS and the recent Super-Kamiokande results from the 
combined phases I, II and III. In Sec. 4 we combine all the data together and we determine 
the presently allowed ranges of mass and mixing, thus updating our previous results. 

2. Leading Am^ oscillations: solar and KamLAND data 

In the analysis of solar neutrino experiments we include the total rates from the radiochemi- 
cal experiments Chlorine [13], Gallex/GNO [14] 1 and SAGE [16]. For real-time experiments 
we include the 44 data points of the electron scattering (ES) Super-Kamiokande phase I 
(SK-I) energy-zenith spectrum [17] and the data from the three phases of SNO [18-20], 
including the results on the low energy threshold analysis of the combined SNO phase I and 
phase II [21] (which we label SNO-LETA). We also include the main set of the 192 days 
of Borexino data [22] (which we label Borexino-LE) as well as their high-energy spectrum 
from 246 live days [23] (Borexino-HE). We consider the following variations on the analysis: 

x Note that the value -RGaiiex = 73.4lg'J (stat) I41 (syst) SNU presented in [14] fully coincide with the 
preliminary result quoted in [15], and has been commonly used in the literature since 2008. 
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Figure 1: Allowed parameter regions (at 90%, 95%, 99% and 99.73% CL for 2 d.o.f.) from 
the combined analysis of solar data for 13 = 0. The best-fit point is marked with a star. For 
comparison we also show as empty regions (the best-fit is marked by a circle) the results prior to 
the inclusion of the latest Ga capture rate of SAGE [16], the energy spectrum of Borexino [22,23] 
and the low energy threshold analysis of the combined SNO phase I and phase II [21]. In both 
analysis we use as inputs the GS98 solar model fluxes and the Gallium capture cross-section of 
Bahcall [24]. 

• Updated capture rate in gallium. In Ref. [16] the Russian- American experiment 
SAGE has presented the results of the combined analysis of 168 extractions (until 
December 2007). The capture rate of solar neutrinos above the 233 keV threshold is 

#sage = 65.4tjj;i (stat) ±i;l (syst) SNU. (2.1) 

which is slightly lower (but fully consistent) than the previous result quoted in 
Ref. [15], 66.2+32 (stat)+3 2 (syst) SNU, and presents a considerable improvement 
of the systematic uncertainties. 

• Possible modification of the capture cross-section in gallium. Ref. [16] also 
presents the results of a new calibration of the SAGE detector with a reactor- produced 

Ar neutrino source. The ratio of observed to expected event rate in this experi- 
ment, once combined with the measured rates in the three prior 51 Cr neutrino-source 
experiments with Gallium, is 0.87 ± 0.05. As a possible explanation for this low 
result, in Ref. [16] it is proposed that the cross-section for neutrino capture by the 
two lowest-lying excited states in 71 Ge may have been overestimated in Ref. [24]. 
As an alternative, the authors consider a modified capture cross-section where the 
contribution from these two excited states is set to zero. 
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• Inclusion of the Borexino spectral data. Following the procedure described in 
Ref. [25] (see Appendix A of that work for details) we have included in the analysis 
the 160 data points of the Borexino energy spectrum in the 365-2000 keV energy 
range [22] as well as the 7 points of the high-energy spectrum from the 246 live days 
of Borexino [23]. In our approach the overall normalizations of the 11 C, 14 C, 210 Bi 
and 85 Kr backgrounds are introduced as free parameters and are fitted to the data. 

• Low energy threshold analysis of the combined SNO phase I and phase II 
data. In Ref. [21] the SNO Collaboration reported the results from a joint analysis 
of their Phase I and Phase II data with an effective electron kinetic energy threshold 
of T e ff = 3.5 MeV. Besides the inclusion of the lower energy data the analysis present 
improvements in the calibration and analysis techniques to reduce the threshold and 
increase the precision of the results. 

From the point of view of any reanalysis of their data, there is an important difference 
with respect to the previous higher threshold results. In Refs. [18,19] the results were 
quoted in the form of binned (in energy and time) event rates. In particular, 34 data 
points of the day-night spectrum were given for SNO-I, and the separate day and 
night rates for neutral current (NC) and ES events as well as the day-night energy- 
spectrum for charge current (CC) events were given for SNO-II. Instead, in Ref. [21] 
the collaboration presents its reanalysis as a total 8 B neutrino flux plus an effective 
description (under the assumption of unitarity for active neutrinos) of the v e survival 
probability, whose dependence on E v is parametrized as a quadratic function for P«f e ay 
and a linear function for the day/night asymmetry. In other words, the SNO-LETA 
results are given as best-fit, uncertainties and correlations of six effective parameters: 
the 8 B flux, three polynomial coefficients of E v for the Pee y and two for the day- 
night asymmetry. When using these data to perform a fit to neutrino oscillations, 
for each value of the oscillation parameters one must first obtain the polynomial 
coefficients that best represent the corresponding Mikheev, Smirnov and Wolfenstein 
(MSW) [26, 27] oscillation probability. This has to be done taking into account 
the sensitivity of the SNO detector, and the relevant information is provided in the 
Appendix A of Ref. [21]. We have followed the procedure outlined there and we have 
verified that we can perfectly reproduce their oscillation results in Fig. 37. 2 

• Uncertainties on the Solar Model. Recent detailed determination of the abun- 
dances of the heavy elements on the solar surface [28, 29] lead to lower values than 
previous studies [30]. Any solar model which incorporates such lower metallicities 
fails at explaining the helioseismological observations [31]. Changes in the Sun mod- 
eling, in particular of the less known convective zone, are not able to account for this 
discrepancy [32,33]. So far there has not been a successful solution of this puzzle, so 

2 However we notice that this procedure does not allow the inclusion of the SNO low energy threshold 
data in analysis in terms of more exotic scenarios in which either unitarity in the active neutrino sector does 
not hold (like for scenarios with sterile neutrinos) or the energy dependence of the oscillation probability 
cannot be well represented by a simple quadratic function. 
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that at present there is no fully consistent Standard Solar Model (SSM). This lead to 
the construction of two different sets of SSM's, one (labeled "GS") based on the older 
solar abundances [30] leading to high metallicity, and one (labeled "AGS" ) assuming 
lower metallicity [34]. We use the most recent recalculation of the fluxes in these 
two sets from Ref. [35], and following their notation we refer to them as GS98 and 
AGSS09. 

Let us first consider the case #13 = 0. In Fig. 1 we show the present determination of the 
leading parameters Am^ and 9i 2 from the updated oscillation analysis of the solar neutrino 
data described above in the context of the GS98 solar model. For comparison we also show 
the results obtained prior to the inclusion of the latest Ga capture rate of SAGE [16], the 
energy spectrum of Borexino [22, 23] and the SNO-LETA results [21] for the same solar 
model. As seen in this figure, the inclusion of these results lead to an improvement on 
the determination of both O12 and Am^ and for this last one the best-fit value slightly 
increases. The effect of the different variations in the analysis is displayed in Fig. 2, where 
we show the 99% CL region in the (Am^, tan 2 612) plane from the present analysis and 
the corresponding one when some of the input data or assumption are modified. The figure 
shows that the most quantitatively relevant new information arises from the inclusion of 
the SNO-LETA results. The inclusion of Borexino tends to shift the region towards slightly 
lower values of 612 angle. Conversely, if the analysis is done in the context of the AGSS09 
model the region is shifted towards slightly larger 612- 

2.1 Impact of #13 on the solar analysis 

The survival probability of solar neutrinos in the framework of three neutrino oscillations 
can be written as: 

P% = sin 4 6 13 + cos 4 9 13 P%(Am 2 21 ,e 12 ) , (2.2) 

where we have used the fact that Lg| c = 47r£ , /Am| 1 is much shorter than the distance 
between the Sun and the Earth, so that the oscillations related to L^\ c are averaged. In 
presence of matter effects P^^Am^, 0± 2 ) should be calculated taking into account the 
evolution in an effective matter density n| = n e cos 2 #13. For 1(T 5 < Am 2 /eV 2 < 1(T 4 , 
Pg^Am^, 612) presents the following asymptotic behaviors [36]: 

P e 2 ;(Am^, 0i 2 ) ~ 1 - * sin 2 (2#i 2 ) for E v < few x 100 KeV (2.3) 

P^(Am 2 21 ,e 12 ) ~ sin 2 (#i 2 ) for E v > few x 1 MeV (2.4) 

The impact of the inclusion of a non-zero value of #13 in the solar analysis is shown 
in Fig. 3 and in the upper-left panel of Fig. 4. From this last plot we see that solar 
neutrino data by themselves favor #13 = 0, although their sensitivity is very weak for 
sin 2 6*i3 < 0.03. This behavior can be understood from Fig. 3, where we show the allowed 
regions (at 95% CL) in the (Am^, tan 2 612) plane as obtained from the analysis of low- 
energy (radiochemical and Borexino-LE) and high-energy (SK, SNO and Borexino-HE) 
solar experiments, for different values of #13. As described in Eq. (2.2), for fixed values of 
Am 2 ] and 612, the inclusion of a small value of #13 results into a decrease on the predicted 
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Figure 2: Effects of the different assumed inputs on the 99% CL parameter region from the 
combined analysis of solar data for #13 = 0. The full region corresponds to the results show in 
Fig. 1. The solid line corresponds to the modification due to the change on one of the inputs as 
labeled in each panel. 

rates at a given solar neutrino experiment. This decrease can be compensated by a shift 
of Am^i an d 012 which lead to an increase of P%£. However the sign of the shift strongly 
depends on the characteristic energy of the detected neutrinos. For experiments detecting 
neutrinos with energies low enough for matter effects to be irrelevant (such as Chlorine and 
Gallium experiments) P^ is given by Eq. (2.3) and increases as 812 decreases. Conversely, 
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Figure 3: Dependence on 9 13 of the allowed (Am^, tan 2 6*12) regions (at 95% CL for 2 d.o.f.) from 
the partial analysis of the "low energy" and "high energy" solar neutrino data. The left column 
corresponds to the analysis prior to the inclusion of the latest Ga capture rate of SAGE [16], the 
energy spectrum of Borexino [22,23] and the low energy threshold analysis of the combined SNO 
phase I and phase II [21]. In the central and right columns the results of those experiments are 
included. The central column corresponds to GS98 solar model fluxes and Gallium capture cross- 
section of Bahcall [24]. In the right panels the AGSS09 solar model fluxes are used. The full (void) 
regions are obtained with Gallium capture cross-section of Bahcall [24] (modified cross-section in 
Ref. [16]). 



for experiments detecting neutrinos mostly in the regime of adiabatic matter oscillations 
(such as SK and SNO) is given by Eq. (2.4) and increases as #12 increases. Consequently 
the combined fit worsens with #13. 

Fig. 3 also illustrates how the inclusion of the new data strengthens this tension. The 
quantitative improvement on the corresponding bound on #13 is displayed in the upper-left 
panel of Fig. 4. We have also studied the dependence of this result on the assumed solar 
model and on the possible modification of the neutrino capture in gallium. We find that 
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Figure 4: Dependence of A% 2 on sin 2 #13 for the different analysis as labeled in the figure and 
the corresponding la ranges. 



the AGSS09 fluxes lead to a stronger bound on #13 and the same happens with the modified 
gallium cross-section. This second effect can be easily understood as follows: because of 
the lower cross-section, a slightly larger survival probability is required to fit the same 
data. This shifts the (Am| l5 tan 2 #12) allowed region obtained from the analysis of the 
radiochemical solar experiments towards slightly lower values of #12 (as explicitly shown 
in the last column of Fig. 3). As a consequence it increases the tension with the SNO+SK 
favored mixing angle for non-zero #13. 
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Figure 5: Allowed parameter regions (at 90%, 95%, 99% and 99.73% CL for 2 d.o.f.) from 
analysis of KamLAND (full regions with best-fit marked by a star) and solar (void regions with 
best-fit marked by a dot) data for two values of #13 as labeled in the figure and for the two solar 
models. 



2.2 Combination with KamLAND and the hint of #13 / 

We show in the left panels of Fig. 5 the present determination of the leading parameters 
Am^ and #12 (for #13 = 0) from the analysis of KamLAND spectral data [37] compared to 
those from the updated solar analysis for the two solar models considered. While the results 
show perfect agreement in Am| 1; there appears to be a mismatch in the favored value of 
012 as determined from KamLAND compared to the one from solar neutrinos, this last one 
being mostly sensitive to the precise value of CC/NC event (i.e., to (P ee ) oc sin 2 6*12) as 
determined by SK and SNO. 

As it was pointed out in Ref. [11] and widely discussed in the literature [8-10,38], this 
mismatch can be lifted by a non-zero value of #13. This happens because, as discussed 
above, the CC/NC event rate can be fitted with a higher value of #12 provided that a 
non-zero #13 is included. Conversely for KamLAND Eq. (2.2) also holds with 

P% Mm = l-\ sin 2 (20 12 ) sm 2 . (2.5) 
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Figure 6: Allowed parameter regions (at 90%, 95%, 99% and 99.73% CL for 2 d.o.f.) from the 
analysis of atmospheric data (full regions, best-fit marked with a star) and LBL data (void regions, 
best-fit marked by a circle) for #13 = and Am|i = 7.6 x 10~ 5 eV 2 . 



So for > the KamLAND spectrum can be well fitted with a smaller value of 612, and 
consequently the best-fit values for solar and KamLAND analysis agree better for #13 ^ 0. 
This behavior is clearly visible in the right panels of Fig. 5 and in the upper-central panel 
of Fig. 4. As mentioned before, the best-fit value of #12 for solar neutrino fit within the 
AGSS09 model is slightly larger than for the GS98 model, and therefore the required value 
of #13 to achieve agreement with KamLAND is smaller for the AGSS09 model. 

However, one must notice that the better agreement between the solar and KamLAND 
analysis for #13 7^ has to be contrasted with the worsening of the global description of 
the solar neutrino data previously described. As a consequence of this tension we find that 
the inclusion of the new data and of the modified gallium capture cross-section also tend 
to lower the best-fit value of #13 as well as its corresponding statistical significance (see 
upper-central panel of Fig. 4). Altogether we find that the la range for #13 as determined 
from the global analysis of solar and KamLAND data changes from the old value sin 2 #13 = 
0.025 ± 0.018 (before the inclusion of SNO-LETA, SAGE-09 and Borexino data) to: 



sin 2 #13 



0.021 ±0.017 for GS98, 

0.017 ±0.017 forAGSS09, (2.6) 
0.015 ± 0.017 for AGSS09 with modified Ga cross-section. 



3. Leading Am^ oscillations: atmospheric, CHOOZ and accelerator data 

In this section we present two different analyses of the atmospheric data. The first one 
is very similar to the one detailed in Ref. [3], and includes the results from the first run 
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of Super-Kamiokande, which accumulated data from May 1996 to July 2001 (1489 day 
exposure) and is usually referred as SK-I [39], as well as the data obtained with the par- 
tial coverage after the 2001 accident (804 day exposure), the so-called SK-II period [40]. 
We will refer to this analysis as SK(I+II). The second one is based on the new analysis 
recently presented by the Super-Kamiokande collaboration including also the data taken 
from December 2005 to June 2007, usually referred as SK-III [41] . Apart from the inclusion 
of these new event rates, in this data release the selection criteria and the corresponding 
estimate of uncertainties for the SK-I and SK-II periods have been changed with respect 
to the previous SK(I+II) analysis. We have therefore performed a reanalysis of the new 
combined samples from phases I, II and III as presented in [41]. We refer to the results 
of this analysis as SK(I+II+III). It is important to point out that already since SK-II the 
Super-Kamiokande collaboration has been presenting its experimental results in terms of a 
large number of data samples. The rates for some of those samples cannot be theoretically 
predicted (and therefore include in a statistical analysis) without a detailed simulation of 
the detector which can only be made by the experimental collaboration itself. Thus our 
results represent the most up-to-date analysis of the atmospheric neutrino data which can 
be performed outside the collaboration. For details on our simulation of the data samples 
and the statistical analysis see the Appendix of Ref . [3] . 

For what concerns LBL accelerator experiments, we combine the results on disap- 
pearance from K2K [42] with those obtained by MINOS at a baseline of 735 km after a 
two-year exposure to the Fermilab NuMI beam, corresponding to a total of 3.36 x 10 20 
protons on target [43]. We also include the recent results on — > v e transitions based on 
an exposure of 7 x 10 20 protons on target [44,45]. 

In order to test the description of the present data in the absence of ^13-induced effects 
we show in Fig. 6 the present determination of the leading parameters Am 2 ^ and #23 for 
#13 = and Anijj = 7.6 x 10 -5 eV 2 from the two atmospheric neutrino analyses and the 
LBL accelerator results. For concreteness we plot only normal ordering; the case of inverted 
ordering gives practically identical results as long as 6*13 = 0. This figure illustrates how 
the bounds on the oscillation parameters #23 and Am^ emerges from a complementarity of 
atmospheric and accelerator neutrino data: [Am^l is determined by the spectral data from 
MINOS, whereas the mixing angle #23 is still largely dominated by atmospheric data from 
Super-Kamiokande with a best-fit point close to maximal mixing. It is important to note 
that there is a very good agreement in the location of the best-fit points from SK(I+II) 
and MINOS. This is not the case for SK(I+II+III) for which the best-fit point in [Am^l 
is now lower than the one obtained from LBL. In what follows we will discuss the impact 
of this tension on the allowed range of #13 . 

3.1 Impact of #13 on the atmospheric and LBL f M disappearance data 

It is well known that a very important contribution to our knowledge of #13 arises from 
the negative results on v e disappearance at short baselines at the CHOOZ reactor experi- 
ment [46] . Given the values of Am^ obtained from solar and KamLAND experiments, the 
disappearance of 9 e at the CHOOZ distances arises from oscillations due to [Am^) whose 
amplitude is proportional to sin 2 {26 13). Consequently, when combined with the bounds on 
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Figure 7: Zenith distribution for multi-GeV e-like events (upper panels), and A\ 2 dependence 
on cosc>cp sin 0i3 marginalized over all undisplayed parameters (lower panels). The fits include the 
CHOOZ data as well as the full atmospheric data samples for SK-I (left panels), SK-I with modified 
multi-GeV e-like data (central-left panels), SK(I+II) (central-right panels) and SK(I+II+III) (right 
panels) . 



jAm^j from atmospheric and disappearance at LBL experiments the CHOOZ result 
implies an upper bound on #13. 

In Ref. [11] a hint for a non-zero value of #13 was obtained at 0.9<t from the analysis of 
atmospheric data from SK-I, z/,j disappearance at long-baseline experiments and CHOOZ. 
The authors traced its origin to subleading effects driven by Am^. This result was ex- 
tensively studied in Ref. [10], where it was pointed out that the appearance of the hint 
is triggered by an excess of multi-GeV e-like data in the first two angular bins of SK-I 
data, which lead to a better fit for this sample in the case of a non-zero value of #13. This 
conclusion is summarized in Fig. 7: as clearly visible in the leftmost panel, there is a small 
but definite preference for a non-zero value of #13 in SK-I. Once the SK-II data are included 
in the analysis the event rates of those bins are considerably reduced, and as a consequence 
no hint of non-zero $13 is found in SK(I+II). As can be seen in the rightmost panels, in 
the new SK(I+II-|-III) data release there is again a small excess in the first angular bin 
of of multi-GeV e-like data, and consequently a weak preference for non-zero #13 reap- 
pears. Such "hint" , however, is considerably weaker than in the SK-I case and completely 
negligible from the statistical point of view. 

In order to better understand the impact of the various data samples, in the second 
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row of Fig. 4 we plot the bound on #13 as implied by different combinations of atmo- 
spheric and LBL experiments. In all panels the results are marginalized with respect to 
the undisplayed parameters, and in particular to the sign of Am^. As can be seen, for 
both SK(I+II) and SK(I+II+III) the combined analysis of atmospheric and CHOOZ show 
no relevant preference for nonzero #13, and the data results exclusively into a bound on #13. 
Such bound is weaker for SK (1+11+111) as a result of the lower values of [Am^-J preferred 
by the atmospheric analysis, and gets slightly strengthened by the inclusion of LBL 
disappearance data. Since the survival probability at MINOS and K2K is practically 
insensitive to such small values of 6*13, this improvement arises as an indirect effect driven 
by the better determination of lAm^j. On the other hand, it should be noted that when 
combining the new SK(I+II+III) with LBL disappearance data a small preference for 
a non-zero value of $13 arises. This effect is a consequence of the mismatch between the 
best-fit lAm^l in atmospheric and MINOS data shown in Fig. 6. Within our atmospheric 
neutrino analysis this is only a A% 2 = 0.5 effect. However, the analysis of the atmospheric 
data presented in Ref. [41] by Super-Kamiokande seems to exhibit a non-negligible depen- 
dence of the preferred Am^ range on #13 (see, e.g., the leftmost panels of Figs. 7 and 8 
in [41]). It would therefore be very interesting if the Super-Kamiokande and MINOS col- 
laborations could perform a combined analysis of their samples in order to fully establish 
the statistical significance of this effect. 

3.2 — > v e appearance results in MINOS 

In Ref. [47] the MINOS collaboration reported their first results on the search for — > u e 
transitions based on an exposure of 3.14 x 10 20 protons-on-target in the Fermilab NuMI 
beam. They observed 35 events in the Far Detector with a background of 27 ± 5 (stat) ± 
2 (syst) events predicted by their measurements in the Near Detector. This result corre- 
sponded to a 1.5<r excess which could be explained by a non-zero value of #13. Recently a 
new analysis with double statistics (exposure of 7 x 10 20 ) has been presented [44,45]. The 
MINOS collaboration reported the observation of 54 events with an expected background 
of 49.1 ± 7.0 (stat) ± 2.7 (syst), thus reducing the excess above background to 0.7a. 

In the left panel of Fig. 8 we show the reconstructed energy distribution of the observed 
events together with the background expectations. We compare our simulation of the 
expected spectrum in the presence of — > v e oscillations with the one given by the 
MINOS collaboration. In the right panel of Fig. 8 we perform an oscillation analysis of the 
MINOS total event rate, and again we compare our limits on #13 for given #cp to those 
obtained by the MINOS collaboration. In both cases we find perfect agreement. 

In the central panel of Fig. 4 we show the information on #13 from the analysis of the 
LBL data from both v e appearance and disappearance and its impact when combined 
with the atmospheric and CHOOZ data. As seen in this figure, the inclusion of MINOS ap- 
pearance data in combination with either atmospheric neutrino samples leads to a stronger 
upper bound on #13. Due to the 0.7<r excess it also increases slightly the significance for 
013 7^ 0. As expected, this increase is larger when combined with SK(I+II+III), since in 
this case the MINOS v e appearance excess adds to the Am^ mismatch effect between 
SK(I+II+III) and MINOS disappearance discussed above. Altogether we find that the 
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Figure 8: Left: reconstructed energy distribution of the v e charged current events in the MINOS 
far detector. We show the background prediction (thick solid green line) as well as our (red solid) and 
MINOS (blue dashed) predictions for the expected spectrum in the presence of — > v e oscillations 
with sin 2 26*13 = 0.115 and Sep = 0. Right: best-fit and upper bound on 6*13 as a function of i5cp for 
the MINOS v e appearance data. The dashed (dotted) lines represent the best-fit as obtained from 
our (MINOS) analysis for normal and inverted ordering. The solid regions are our 90% bounds 
for normal (darker blue) and inverted (lighter red) ordering; the solid lines show the corresponding 
limits from MINOS, taken from Ref. [44]. In both panels the remaining parameters are fixed to 
tan 2 6» 12 = 0.45, maximal <9 23 , Am^ = 7.6 x 1(T 5 eV 2 and |Am 2 2 | = 2.43 x 10~ 3 eV 2 . 

la range for #13 as determined from the analysis of atmospheric, CHOOZ and LBL data 
reads: 

sin 2 = O.0O2±g;g5| CHOOZ + LBL(DIS+APP) + SK(I+II) , (3.1) 
sin 2 0i3 = 0.00618:8$ CHOOZ + LBL(DIS+APP) + SK(I+II+III) , (3.2) 

where for the sake of illustration we have extrapolated the lower la bound to the unphysical 
region sin 2 #13 < 0. 

4. Global results and conclusions 

The results of the global combined analysis including the SK(I+II+III) atmospheric neu- 
trino data and all dominant and subdominant oscillation effects are summarized in Fig. 9, 
where we show the different projections of the allowed six-dimensional parameter space. 
The full regions correspond to the analysis done in the framework of the GS98 solar model 
and with Ga capture cross-section in Ref. [24] while the void regions correspond the anal- 
ysis with AGSS09 solar fluxes and the modified Ga capture cross-section in Ref. [16]. The 
regions in each panel are obtained after marginalization of Xgiobal with respect to the undis- 
played parameters. In the lower panel we show the allowed regions in the (sin 2 $13, <5cp) 
plane. As seen in the figure, at present the sensitivity to the CP phase is marginal but we 
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Figure 9: Global Sv oscillation analysis. Each panels shows two-dimensional projection of the 
allowed five-dimensional region after marginalization with respect to the undisplayed parameters. 
The different contours correspond to the two-dimensional allowed regions at 90%, 95%, 99% and 
3a CL. The full regions correspond to the analysis done in the framework of the GS98 solar model 
and with Ga capture cross-section in Ref. [24] while the void regions correspond the analysis with 
AGSS09 solar fluxes and the modified Ga capture cross-section in Ref. [16]. 
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find that the bound on sin 2 #13 can vary by about ~ 30% depending on the exact value of 
5cp. This arises mainly from the interference of #13 and Am^ effects in the atmospheric 
neutrino observables, as well as from the new MINOS v e appearance data. The derived 
ranges for the six parameters at the la (3d) level are: 

GS98 with Gallium AGSS09 with modified 

cross-section from [24] Gallium cross-section [16] 

Amli = 7.59 ± 0.20 (±g;jj£) x 10" 5 eV 2 Same 

„ f— 2.36 ± 0.11 (±0.37) x 10~ 3 eV 2 
Ami-, = < 00 Same 

[+2.46 ± 0.12 (±0.37) x 10~ 3 eV 2 

012 = 34.4 ± 1.0 (1| 2 )° 34.5 ± 1.0 (^ 3 ; 2 )° 
6 23 = 42.8 tt 7 9 ( +W 7 l)° Same 

13 = 5.6 t 3 ;° (< 12.5)° 5.1 ± 3 3 ° 3 (< 12.0)° 

[ sin 2 #13 = 0.0095 (< 0.047) ] [0.008 (< 0.043) ] 

(5 C p G [0, 360] Same 

For each parameter the ranges are obtained after marginalizing with respect to to the 
other parameters. For Am^ the allowed ranges are formed by two disconnected intervals 
which correspond to the two possible mass orderings. The absolute best-fit lays in the 
positive Am 2 ! = ±2.46 x 10~ 3 eV 2 . The la and 3a ranges are defined with respect 
to this absolute minimum. In particular the local best-fit in the inverse mass ordering, 
Am§! = -2.36 x 10~ 3 eV 2 , is at a Ax 2 = 0.12. 

In summary we have presented the results of an up-to-date global analysis of neutrino 
oscillation data including the detailed analysis of the Borexino spectra, the lower value of 
Ga rate as lastly measured in SAGE, the low energy threshold analysis of the combined 
SNO phase I and phase II, the results for v e appearance from MINOS and the recent 
reanalysis of SK(I±II±III) atmospheric data. We have studied the robustness of the hints 
of a non- vanishing value of #13 under the inclusion of the new data and variations of the 
assumptions about the solar model and the possible modification of the cross-section for 
neutrino capture in Ga as advocated by the SAGE collaboration in order to explain their 
new calibration data. 

We found that the inclusion of the new solar data, and in particular of the SNO- 
LETA results tends to lower the statistical significance of #13 7^ while the results from 
v e appearance from MINOS and the new SK(I±II±III) atmospheric neutrino reanalysis 
increases it. We conclude that the significance of #13 7^ is: 

GS98 with Gallium AGSS09 with modified 

cross-section from [24] Gallium cross-section [16] 

solar±KamLAND CL = 79% (1.26<r) CL = 70% (1.05a) 

x ^L(DiJSi) CH00Z cl = 8i% (l3ict) cl = 76% (Li7a) 

The determination of the other oscillation parameters is rather robust. 
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